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Abstract —We consider the use of multiple co-located satellites 
to improve the spectral efficiency of broadcast transmissions. In 
particular, we assume that two satellites transmit on overlapping 
geographical coverage areas, with overlapping frequencies. We 
first describe the theoretical framework based on network in¬ 
formation theory and, in particular, on the theory for multiple 
access channels. The application to different scenarios will be 
then considered. Including the bandlimited additive white Gaus¬ 
sian noise channel with average power constraint and different 
models for the nonlinear satellite channel. The comparison with 
the adoption of frequency division multiplexing and with the 
Alamouti space-time block coding is also provided. The main 
conclusion is that a strategy based on overlapped signals is the 
most convenient in the case of no power unbalance, although it 
requires the adoption of a multiuser detection strategy at the 
receiver. 

Index Terms —Co-located satellites. Frequency division multi¬ 
plexing, Multiple access channel. Spectral efficiency. 


I. Introduction 

In today’s satellite communication systems, the scarcity of 
frequency spectrum and the ever growing demand for data 
throughput has increased the need for resource sharing. In 
recent years, users of professional broadcast applications such 
as content contribution, distribution, and professional data 
services have demanded more spectrally efficient solutions. 

Satellite service providers often have the availability of 
co-located satellites: two (or more) satellites are said to be 
co-located when, from a receiver on Earth, they appear to 
occupy the same orbital position. Co-location of satellites 
is typically used to cover the fully available spectrum by 
activating transponders on different satellites that cover non¬ 
overlapping frequencies or as a stand-by in-orbit redundancy, 
when the backup satellite is activated in case of failure of 
the main satellite. However, the second satellite can also be 
exploited to try to increase the capacity of the communications 
link. 

In this paper, we address a scenario in which the backup 
satellite is activated in addition to the main one, to improve the 
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spectral efficiency (SE) of the overall communication system. 
The transmission from the two satellites can be coordinated, 
but through simple geometrical considerations it can be easily 
shown that even with a coverage area of a few tens of 
kilometers and two co-located signals separated in angle by 
a fraction of degree, time alignment is not possible. On the 
other hand, the considered system model can also represent 
a scenario where a single satellite with two transponders 
operating at the same frequency is employed and hence the 
two transmitted signals can be considered synchronous. 

Here, we study the information rate (IR) achievable by 
a system where the two satellites transmit on overlapping 
geographical coverage areas with overlapping frequencies, and 
compare our results with that achievable by the frequency 
division multiplexing (EDM) strategy and with that achievable 
adopting the well-known Alamouti space-time block code |[I]. 

The two-satellites scenario has been studied in H), 0, 
where the satellite channel is approximated as a linear additive 
white Gaussian noise (AWGN) channel, and the information 
theoretic analysis has been carried out under the limiting 
assumption of Gaussian inputs. We instead examine three 
different models for the system: the linear AWGN channel, the 
peak-power-limited AWGN channel 0-0, and the satellite 
channel adopted in the 2nd-generation satellite digital video 
broadcasting (DVB-S2) standard 0. The studied system is 
an instance of broadcast channel [sLfTO) with multiple trans¬ 
mitters. However, we are interested in a scenario in which 
the same information must be sent to every receiver. This 
situation corresponds, for example, to the delivery of a TV 
broadcast channel. We show that all these scenarios can be 
analyzed by means of network information theory and we 
will resort to multiple access channels (MACs) 0^113 with 
proper constraints. 

Our analysis reveals that, if we allow multiuser detection, 
the strategy based on overlapping signals achieves higher SEs 
with respect to that achievable by using EDM. Interestingly, 
we show that there are cases in which a single satellite can 
outperform both these multiple satellites strategies, but not the 
Alamouti scheme. 


The remainder of this paper is structured as follows: in 
Section we present a general system model valid for all 
cases, and in Section we briefly review the theory of 
MACs. In Sections IIVI and lYl we discuss the achievable 
rates by EDM and by the Alamouti space-time code. In 
and 


Sections VI VII 


VIII we analyze the three different 


scenarios and Section IX concludes the paper. 
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II. System Model 

Fig. depicts a schematic view of the baseband model we 
are considering. A single operator properly sends two separate 
data streams to the two satellites. The impact of the feeder 
uplink interference is considered negligible in this scenario. 
Data streams are linearly modulated signals, expressed as 

= - kT) 2 = 1 , 2 , ( 1 ) 

k 

where is the fc-th symbol transmitted on data stream i, 
p{t) is the shaping pulse, and T is the symbol time. 

Each satellite then relays the signal, denoted as Sz(f)QtO 
several users scattered in its coverage area. For each user, the 
received signal is the sum of the two signals coming from the 
satellites, with a possible power unbalance 7 ^ due to different 
path attenuations (we assume 1/2 < 7 < 1). Without loss 
of generality, we assume that the attenuated signal is S2{t), 
otherwise we can exchange the roles of the two satellites. The 
received signal is also affected by a complex AWGN process 
w{t) with power spectral density Nq. As mentioned, time 
alignment between the signals transmitted by the two satellites 
is not possible, since if the signals from the two satellites come 
perfectly aligned at a given receiver in the area, there will be 
other receivers for which a misalignment of a few symbols is 
observed. On the other hand, it is straightforward to show that 
our information-theoretic analysis does not depend on the time 
alignment of the two signals, and we will assume synchronous 
users to simplify the exposition. Hence, the received signal has 
the following expression 

y{t) = si{t) + + w{t) , ( 2 ) 

where Si{t) and S2{t) are the signals at the output of the 
two satellites, and (j){t) is a possible phase noise process, 
caused by the instabilities of the oscillators. We assume that 
the phase noise is slow varying w.r.t. the signals’ baud rate 
and perfectly known at the receiver. Signals si{t) and S2{t) 
are transmitted with overlapping frequencies, and the overall 
signal has bandwidth W. 

Since we are analyzing a broadcast scenario in which 
different receivers experience different (and unknown) levels 
of power unbalance, we impose that the two satellites transmit 
with the same rate. This constraint will be better clarihed in 
the next sections. Channel state information is not available at 
the transmitter and no cooperation among the users is allowed. 
This is because the target is on broadcasting applications. 

A simple alternative strategy to overlapping frequencies, 
that allows to avoid interference between the two transmitted 
signals, is EDM. The bandwidth W is divided into two 
equal subbands assigned to the different satellites. An unequal 
subband allocation does not make sense since the power 
unbalance is different for different receivers in the coverage 
area and, in any case, unknown to the transmitter. In this case, 
the received signal has expression 

y{t) = + w{t) , (3) 

’it can be Si{t) A ^i{i) due to the nonlinear transformation at the satellite 
transponder. 



Figure 1. Block diagram of the analyzed system. 


where fc is the frequency separation between the two signals. 

Another possible alternative to avoid interference between 
the two signals is the use of the Alamouti space-time block 
code Q, consisting in the two satellites exchanging the trans¬ 
mitted signals in two consecutive transmissions. Unlike the 
two previous strategies, its classical implementation requires 
a perfect alignment in time of the signals received from the 
two satellites. However, in Appendix we will describe an 
alternative implementation working in the presence of a delay 
which can be different for different receivers. 

In the following, these transmission strategies will be com¬ 
pared by using the overall SE of the system as a hgure of 
merit. The SE is dehned as 

SE=^ [bit/s/Hz], 

where I is the maximum mutual information of the channel. 
However, since in the scenario of interest the values of T 
and W are fixed, without loss of generality we will assume 
that TW — 1 and we will refer to the terms IR and SE 
interchangeably. 

III. Multiple Access Channels 

In this section, we review some results on classical 
MACs Q. We consider the transmission of independent 
signals from the two satellites|^ We denote by i?i the SE of 
the hrst satellite and by R 2 that of the second satellite. At this 
point, we make no assumptions on the channel inputs, since a 
better characterization of the input distributions is presented in 
the next sections. However, independently of the form assumed 
by the input distribution, the boundaries of the SE region can 
be expressed, for each fixed signal-to-noise ratio (SNR), as Q 

Ri < I{xi-y\x2) = h 
R2 < I{x 2 \y\xi) = I2 
Rl+ R2< I{Xi,X2] y) = /j , 

where I{xi\y\x 2 ), I[x 2 ]y\xi) and I{xi,X 2 \y) represent, re¬ 
spectively, the mutual information between xi and y condi¬ 
tioned to a; 2 , that between X 2 and y conditioned to Xi and 
that between the couple (a;i,a: 2 ) and y\ we have omitted the 
dependence on t and we adopt dehnitions Ji, I 2 , and /j to 
simplify the notation. 

^We will explain later in Section |Vl| why this is the best choice for the 
signals transmitted from the two satellites. 
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Figure 2. Achievable rate region in the case I 2 < I\- Point D is the 
maximum achievable SE from satellite 1 to the receiver when satellite 2 is 
not sending information. Point C is the maximum rate at which satellite 2 can 
transmit as long as satellite 1 transmits at its maximum rate. The maximum 
of the sum of the SEs is obtained on points of the segment B-C, which can 
be achieved by joint decoding. 

Fig. I^is useful to gain a better understanding of the behavior 
of SE regions. Point D corresponds to the maximum achievable 
SE from satellite 1 to the receiver when satellite 2 is not 
sending any information. Point C corresponds to the maximum 
rate at which satellite 2 can transmit as long as satellite 1 
transmits at its maximum rate f] The maximum of the sum of 
the SEs, however, is obtained on points of the segment B- 
C; these points can be achieved by joint decoding of both 
signals. It does not make sense to adopt different rates for the 
two satellites, since each satellite ignores whether its signal 
will be attenuated or not and this attenuation will vary for 
different receivers. As a consequence, the only boundary point 
of the SE region we can achieve is point E, which lies on the 
line i?i = i?2. We define /j p the pragmatic sum of rates 
corresponding to point E: it is easy to see, through graphical 
considerations, that 

/j,p = min(/j,2/2) . 

Point E is the intersection between the capacity region and the 
line R2 = —+ Ii and it corresponds to a sum-rate equal 
to Ii. The position of point E depends on the power and on 
the power unbalance. Depending on these two values, E can 
be found in different positions: in particular, if E lies on the 
left of E, we can notice that /j p < Ii, hence it is convenient 
to use a single satellite with rate Ii rather than activating the 
second satellite. 

IV. Achievable Rates by Erequency Division 
Multiplexing 

Since the two signals transmitted by the EDM model ([^ 
operate on disjoint bandwidths, they are independent and the 
IR achievable by this system is equal, in case 7 = 1 , to that 
of a single transmitter with double SNR. We define by /fdm 
the achievable rate by EDM, and by /fdm,p that achievable 
by EDM under the equal rate constraint. The latter is clearly 
equal to twice the minimum SE of the two subchannels. We 

^If we exchange the role of the two satellites, the same considerations hold 
for points A and B instead of D and C. 


demonstrate that the rate achievable by EDM is always lower 
or equal than that achievable with two signals with overlapping 
frequencies in the absence of nonlinear distortions; the same 
result holds for the pragmatic rates and can be stated in the 
following theorem, whose proof can be found in Appendix 

Theorem 1. Let us consider the ideal multiple access channel 

y(t) = xi{t) + ^X2{t) + w(t). ( 4 ) 

The following inequalities hold 

I3 > IpDM ( 5 ) 

.^J,p > ^fdm.p (6) 

with equality if and only if Xi{t) and X2{t) are Gaussian 
random processes and 7^ = 0 dB. 

The theorem, beyond the mathematical proof, has a practical 
explanation. The use of a second satellite, besides increasing 
the overall transmitted power, makes the distribution of Xi {t) + 
7x2 (t) closer to a Gaussian distribution (see the Berry-Esseen 
theorem Thus, a sort of shaping gain must be added to 

the gain arising from the higher power. 

The strategy based on EDM is perfectly equivalent, in 
terms of SE, to a strategy based on time-division multiplexing 
(TDM), in which time is divided into slots of equal length, 
and each satellite is allotted a slot during which only that 
satellite transmits and the other remains silent. During its slot, 
each satellite is allowed to use twice the power. However, on 
satellites, due to peak power constraints, it is not possible to 
double the power and the satellite amplifiers are not conceived 
for power bursts. Hence TDM strategy will not be considered. 

V. Achievable Rates by the Alamouti Scheme 

We now consider the application of the Alamouti 
scheme Q. The two satellites first transmit xi{t) and 0:2(f) 
and then x^it) and —xl{t), respectively. The rate /a, achiev¬ 
able by the Alamouti scheme, satisfies the following theorem, 
proved in Appendix [C| 

Theorem 2. Let us consider the ideal multiple access channel 

y{t) = xi{t) -b yX2{t) -b w{f) , 

where Xi(t), i = 1,2 are random processes such that —Xi{t) 
has the same finite-dimensional distributions as Xi(f). The 
following inequality holds 

Ij > Ia > .i^FDM 

(a) (b) 

with equality in (a) if and only if xi(t) and X2(t) are inde¬ 
pendent Gaussian random processes with the same variance, 
and in (b) if and only if =0 dB. 

Theorem | 2 ] shows that the Alamouti scheme IR is between 
the ones achievable by two overlapping signals and by EDM. 
However, it has the interesting feature that it is not degraded 
by the equal rates constraint, being Ia,p = Ia, where the 
subscript p stands for pragmatic. This is due to the fact that 
both signals are transmitted once by the satellite with no 
attenuation and once by the satellite with attenuation 7. Hence, 
while it is always true that Ia > /fdm,p. it can happen that 
Ia > Ij,p- 
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VI. Additive White Gaussian Noise Channel with 
Average Power Constraint 
A first case study, useful to draw some preliminary con¬ 
siderations about the theoretical limits for the system under 
consideration, is the classical AWGN channel with average 
power constraint. For this case, we have that the two satellites 
of Fig.[T]have no effect on the signal, hence the received signal 
reads 

y{t) = xi{t) + + w{t). 


We express the average power constraint as 

E[|a:,(<)p] <P 1 = 1 , 2 , 

where P is the maximum allowed average power. 

For this channel, the capacity is reached with independent 
Gaussian inputs, p{t) = sinc(f/T), and TW = 1 p^ . A 
sufficient statistic is derived by sampling the output of a low 


pass filter | 13 |, | 14 |. Since we are assuming slow-varying 
phase noise, the observable is 


( 1 ) 

yk = xi' 


■ 'ye 




■ Wk , 


where (j)k = 4 >{kT). The phase noise does not change the 
statistics, and hence the SE Jj is given by the classical 
Shannon capacity, taking into account the total transmitted 
power, and reads 

/j =log2 ^l-b(l-b7^)^^ , 

where N = NqW is the noise power in the considered 
bandwidth. If, instead, we adopt the FDM model (|^, the SE 
can be computed as the average of the SEs of two subchannels, 
each transmitting on half the bandwidth; 

/fdm = ^ log 2 (^1 + 2^^ -b ^ log 2 (^1 + 27^^^ • 

When we introduce the equal rate constraint, it is straightfor¬ 
ward to show that we have the following pragmatic SEs 


In case of EDM, we clearly see how the user with the lower 
SE limits /fdm.p- The curves coincide for 7=1, while they 
suffer from a significant performance loss w.r.t. /fdm for high 
values of power unbalance. If the power unbalance is very 
high, EDM performs even worse than a single satellite. Einally 
we can notice that, when 7^ 7^ 0 dB, /a > /j,p for low P/N 
values. 

At the end of this section, we would like to motivate 
our choice of transmitting independent signals from the two 
satellites. Let us consider the opposite scenario where the 
same signal is transmitted from the two satellites. The received 
signal y{t) can thus be expressed as 

y{t) = x{t) + 'ye^'^^*'^x{t — r) -b w{t) 

where x{t) is the transmitted signal, and t the difference 
between the propagation delays of the two satellites. The 
received sample at time kT reads 

yk = Xk+ 'yel'^’’ ^ sinc(f - T/T)xk-i + Wk ■ 

The channel is thus equivalent to a time-varying frequency- 
selective channel 


Vk 


— ^ ^ kkyXk—i 


■ Wk 


( 7 ) 


with impulse response hk^i = 7ef‘^'=sinc(f — r/T) for f 7^ 0 
and hkfi = 1 + 7ef‘^''sinc(—r/T). As already said, (j)k is 
assumed slowly varying with respect to the symbol interval 
but, due to the oscillators’ instabilities, it will be assumed with 
a coherence time shorter than the codeword length. Hence, we 
are interested in the ergodic rate obtained by averaging the 
information rate that can be obtained with a given value of (j). 
Independently of the value of r, the average signal power is 


E 


^ ^ kkyXk—i 


= (1 + 7^)P 


( 8 ) 


Jj,p = min ^/j,21og2 (^ 1 + 7^^^^ 

7 fdm,p = log2 + 27^ —^ . 

In Eig. we show the SE /j as a function of P/N, 
for different values of power unbalance 7, together with the 
SE that can be achieved when a single satellite is available 
(7 —> 0 ). In this case, the performance of the Alamouti scheme 
is exactly the same as /j, as foreseen by Theorem]^ The figure 
also shows /fdm for the same values of 7. We see that EDM is 
capacity-achieving when 7 = 1 (i.e., /j = /fdm when 7 = 1, 
as also clear from the equations and as foreseen by Theorem[T]) 
but it is suboptimal in the case of power unbalance. 

In Eig. 1 ^ we report the pragmatic SEs for the cases of 
Eig. ID Eor signals with overlapping frequencies, with power 
unbalance 7^ 7/ 0 dB, /j p is lower than /j only in the range 
of low P/N values, corresponding to the case 2/2 < /j. The 
transition is indicated by the change of slope in the curve. We 
also see that, for high power unbalance, a portion of /j p lies 
below single-satellite SE. 


and it can be show that the ergodic rate cannot be higher than 
log2 (l + (l + 7^);^), the rate achievable with independent 
signals (see Appendix [P] for a detailed proof). 

VII. Additive White Gaussian Noise Channel with 
Peak Power Constraint 

As a first step to the theoretical characterization of our 
satellite transmission problem, we consider the case a peak- 
power-limited signal rather than an average-power-limited one. 
The adoption of a peak power constraint comes naturally 
from the use of a saturated nonlinear high-power amplifier 
(HPA) at the satellite. However, there is no expression for 
the channel capacity in this scenario, but only bounds are 
available Q. Eor this reason, we concentrate on the study 
of a simplified discrete-time channel, where the peak power 
constraint is imposed on information symbols ||^. 

In this section, we repeat the analysis of Section |Vl] in 
a peak-power-limited scenario. We first review the results 
in 0 for the case of a single transmitter, then we extend 
the reasoning to the case of two transmitters. 
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Figure 3. Joint spectral efficiency for different values of 7 (AWGN channel 
with average power constraint). 



Figure 4. Pragmatic spectral efficiency for different values of 7 (AWGN 
channel with average power constraint). 


A. Analysis for Single Transmitter 

If we assume that 7 —> 0, the model (|^ simplifies to the 
following discrete-time memoryless channel model 

yk=Xk+wu, (9) 

where is the observable, is the fc-th symbol 

transmitted by satellite 1, and Wk is AWGN with variance 
N = NqW. The input symbols Xk must be subject to a peak- 
power constraint, that can be expressed in the form 

\xk?<P- ( 10 ) 

Channel under constraint ( [TOl l was completely studied 
in 1 ^: the capacity-achieving distribution is discrete in ampli¬ 
tude and uniform in phase, and has the following expression 

^ m 

p{r, 9) = p{9)p{r) = — ^ g^(5(r - ), (11) 

^ e=i 

with Xk = re^^. The distribution is formed of m concentric 
circles, each having weight qg and radius pi. The constraints 


of the problem, in polar coordinates, become 


0 < < y/p 

( 12 ) 

Pi+i > Pe 

(13) 

0 <qe<l 

m 

(14) 

® 1 ■ 

i=i 

(15) 


For the distribution o, we can compute the rate I{xk; Pk) 
in closed form. First of all, we need to derive an expression 
for the probability density functions (PDFs) of the channel and 
the observable. Based on the channel model (|^, we have 

p{yk\xk) = p{yk\r,e) =—e ^ . (16) 

From ([T 6 ]| we can obtain the PDF p{yk) as 


p{yk) 


p+00 


^27r 


p{yk\r, 9)p{r, 9)drd9 


/r=o Je=o 


1 1 r+^ 

ttN 27r 



N 


e ^ 


-pi)drd9 

e=i 


1 1 
ttN 2tt 




1=1 


|yfcP + |p^l^ 

N 


(.271- 


0=0 


^\vk\P( 


cos(arg(yfc)-0)^^ 


1 




£=1 


« lo 


/ 2 |t/fcbA 

1 N ) 


(17) 


where Io( ) is the modified Bessel function of the first kind 
and order zero. By combining ( [T 6 | ) and © we have 


l{xk\yk) = E 


= E 


log; 


p{yk\xk) 

p{yk) 


(18) 


logs 


f^\vdpA 

1^1=1 m ~ lo AT J 


The expectation in ( [TSl l is taken with respect to the actual 
random variables, i.e., Xk (and hence r and 9) and i/k- The 
latter is a function of Xk and Wk, and thus it depends on 
their statistics. The optimal values of m, qi and pi cannot be 
found in closed form, but they are subject to optimization |^. 
For this reason, we evaluated for increasing values of m 
and, for each value, we optimized the m radii to achieve the 
highest IR. Optimization results for 1 < to < 20 are plotted 
in Fig. 1^ We see that, as expected, as P/N increases, the 
optimal distribution is formed of a higher number of circles. 
We also point out that each curve in Fig. is the envelope of 
all curves with a lower number of circles, so to must be read 
as the maximum number of circles, i.e., one or more circles 
can have zero probability. The optimal number of circles is 
shown in Fig. as a function of P/N. We point out that the 
results in Fig.|^differ from those in Q because of a different 
SNR definition. In fact, in Q, capacity curves are computed 
as a function of the SNR per dimension, while our curves are 
a function of the total SNR. 
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Figure 5. Spectral efficiency for single transmitter for distributions composed Figure 7. Joint spectral efficiency for different values of 7 (AWGN channel 

of an increasing number of circles m (AWGN channel with peak power with peak power constraint). 

constraint). 



Figure 6 . Optimal number of circles of the capacity-achieving distribution 
for single transmitter. 


B. Analysis for Two Transmitters 

Aim of this section is to extend the results of Section IVII-AI 
to the case of two transmitters. For this scenario, we make the 
assumption that the optimal distributions of the two inputs are 
still in the form ( [TT] ). This result has been demonstrated for real 
inputs |T5), but not for complex inputs, as the case of interest 
here. For this reason, the computed IR is a lower bound to 
the actual channel capacity, whose expression is not known. 
Under this assumption, the input amplitude distributions are 

rrii 

Pi^^) ^ = 1,2 

e=i 

and the received signal is an extension of SQ 

( 1 ) , ( 2 ) , 

Uk = xl' + 'yxl + Wk 

= + Wk , 

point out that a phase noise term should be considered in the second 
signal. However, since this shift is assumed to be perfectly known at the 
receiver and the input distributions are invariant w.r.t. a phase rotation, we do 
not add it to our model. 


with each of the two inputs satisfying constraints (|T^-([T5]l. 
In this scenario, we can express the joint IR /j = I{xi,X2; y) 
as an extension of (fTS]): 


= E 


logs 




l 2 


2^i=i 2^i=i 


(19) 


where 


27r /‘27r 


Ar,i = 


p 

Jo Jo 






■d6'id^2 • 


/o ^0 

We have computed for different levels of power unbal¬ 
ance between the two received signals; we have verified that 
the best performance is achieved when using input distribu¬ 
tions with only one circle (i.e., mi = m 2 = 1 in ([T^). The 
joint IR is shown in Fig. |7] where the curve labeled 1 satellite 
is obtained as the envelope of the curves of Fig. We report 
here, for comparison, the IR computed when FDM is used, 
assigning half of the bandwidth to each of the satellites. We 
see that, unlike the case of average power constraint, FDM 
is not the optimal choice, not even in the absence of power 
unbalance (when FDM gains exactly 3 dB from the single 
satellite). This result comes from a straightforward application 
of Theorem In effect, since the two input distributions are 
not Gaussian, the inequality Q is strict. Fig. also reports 
the IR /a, achievable by the Alamouti scheme. As foreseen by 
Theorem]^ we see that for 7 ^ = 0 dB the rate I a is perfectly 
equivalent to IfdMj while for 7 ^ = —6 dB FDM performs 
worse. For all values of 7 we have that /j > /a, since the 
input signals are not Gaussian processes. 

We also point out that, unlike what happens when the 
constraint is on the average power, in this case the theoretical 
upper bound for /j, when 7 = 1, is 6 dB higher than the 
single transmitter case. This is because, if the two signals are 
perfectly in phase, the overall signal has double amplitude, 
and hence its power is 4P (i.e., 6 dB higher). This situation 
is unrealistic (and, in fact, we do not experience a 6 dB gain), 
but it is the upper limit for the IR. 

As already mentioned, in a broadcast scenario we have the 
further constraint that the two transmitters must use the same 
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Ri / (bit/s/Hz) 


Figure 8 . Pragmatic spectral efficiency for different values of 7 (AWGN Figure 9. Spectral efficiency regions for 7 ^ = -6 dB. 
channel with peak power constraint). 


rate. When we impose this constraint to the rates shown in 
Fig. 0 we obtain the pragmatic rates in Fig. We see again 
that, as expected, the rates /j p and /fdm.p have suffered a 
degradation for 7 ^ = —6 dB, and we also notice that, for 
low values of P/N and a high power unbalance, the use 
of a single satellite may be convenient over the use of two 
overlapped signals. However, since the Alamouti scheme is 
not degraded by the application of the equal rates constraint, 
we can conclude that the /a grants the best performance in a 
certain range of P/N. 

We can better understand the impact of the equal rates 
constraint on the joint IR by studying the SE regions of 
the channel for different values of P/N, reported in Fig. 
From the analysis of these figures, we can conclude that the 
maximum sum-rate cannot always be achieved and we can 
have a numerical insight of the values of P/N and 7 ^ that 
allow to improve the rates with respect to a case with only 
one transmitter. In particular, it is easy to understand that 
when the power unbalance is low ( 7 ^ —0 dB) the spectral 
efficiency region is perfectly symmetric and maximum sum- 
rate is achieved for all values of P/N. On the other hand, with 
high power unbalance ( 7 ^ = — 6 dB) it is clear that maximum 
sum-rate can be achieved only at high P/N, whereas when 
the power is low the performance of two satellites is worse 
than that of a single satellite. 


C. Practical Constellations for AWGN Channel with Peak 
Power Constraint 


We are now interested in evaluating the performance of 
practical constellations with a finite number of points on 
the AWGN channel with peak power constraint, in order to 
find which kind of discrete constellations can be successfully 
adopted on the satellite channel. 


Starting from the single transmitter case, we see in Fig. 10 


that M-ary phase-shift keying and amplitude-phase-shift key¬ 
ing (MPSK/MAPSK) constellations, usually adopted in satel¬ 
lite communications, are practically capacity-achieving. How¬ 
ever, as foreseen also by the theoretical analysis, constellations 



P/N / (dB) 


Figure 10. Single transmitter capacity and spectral efficiency for PSK/APSK 
constellations. 


with multiple circles (such as APSK) are suboptimal when 


two transmitters are adopted. This can be seen in Fig. 11 


showing the envelopes of the SEs achievable with quaternary 
PSK (QPSK), 8 PSK, 16APSK, 32APSK and 64APSK, where 
APSK exhibits a loss with respect to the bound /j for high 
P/N values. As suggested by the theoretical results, we see 
that the bound is achieved by replacing APSK constellations 
with PSKs with the same cardinality, whose envelopes are 
again shown in the figure. Pig. [T2] reports the same analysis 
for the pragmatic rates, and the same conclusions hold. We 
point out that EDM and Alamouti schemes perform single- 
user operations, so they practically achieve their corresponding 
theoretical bounds with classical PSK/APSK constellations. 
Pinally, we mention that we have attempted an optimization 
of the constellations, using the same algorithm described 
in |T§, imposing that the constellations adopted by the two 
transmitters are identical. Optimization results suggest that 
PSKs are practically optimal in this case. 




































































Figure 13. 


Block diagram of the considered satellite transponder, 
transponder bandwidth 



-/c/2 


/c/2 


signal 

bandwidth 


Figure 11. Joint spectral efficiency for PSK/APSK constellations and Figure 14. Transponder bandwidth allocation for FDM. 
different values of 7 . 



Figure 12. Pragmatic spectral efficiency for PSK/APSK constellations and 
different values of 7 . 


VIII. Satellite Channel 


This section investigates the performance of the system 
in Fig. [T] when a realistic satellite transponder model is 
used. The block diagram of the adopted transponder model 


is depicted in Fig. 13 It shows an input multiplexing (IMUX) 


filter which removes the adjacent channels, a HPA, and an 
output multiplexing (OMUX) filter aimed at reducing the spec¬ 
tral broadening caused by the nonlinear amplifier. The HPA 
AM/AM and AM/PM characteristics and the IMUX/OMUX 
impulse responses are those provided in Q. Particularly, the 
OMUX filter has -3 dB bandwidth equal to 38 MHz and 
the HPA has AM/AM and AM/PM characteristics that are 
called “conventional” in Q. Although the HPA is a nonlinear 
memoryless device, the overall system has memory due to the 
presence of IMUX and OMUX filters. 

The transmitted signals at the input of the two satellites are 
linearly modulated as in Q, with the same pulse and symbol 
interval, and the information symbols are drawn from a 
discrete constellation. The symbol intervals of the two signals 
are also assumed to be perfectly aligned. Thus, the received 


signal reads as in 0. Process </)(/) models the difference 
of phase between oscillators and their phase noise, and is 
considered perfectly known at the receiver. We employ the 
adaptive receiver proposed in jib) , flT) : a sufficient statistic 
for detection is extracted by using oversampling at the output 
of a low pass filter 0’ and a fractionally-spaced minimum 
mean square error (FS-MMSE) equalizer, working at twice 
the symbol rate, acts as adaptive filter followed by a multiuser 
detector. The adaptivity is accomplished by means of the least 
mean square or the recursive least square algorithms p8) . The 
multiuser detector computes the a posteriori probabilities of 
the symbols as 




cx exp 


Vk- 13 


.( 1 ) 




Nn 


where yk is the sample at the output of the FS-MMSE equal¬ 
izer, /3 is a possible (complex-valued) bias, and (pk = 4>{kT) is 
the phase noise process at the receiver (under the assumption 
that (pit) is slow enough w.r.t. the symbol time). 

Similarly to previous sections, we also consider an EDM 
scenario; the transponder bandwidth is equally divided into 
two subchannels as schematically depicted in Eig. 14 Then, 
the EDM receiver performs detection separately with two ES- 
MMSE equalizers, followed by a symbol-by-symbol receiver. 

As already done with the other channel models, we adopt 
the Alamouti scheme as a third possibility: the Alamouti 
precoding is performed on transmitted symbols and, at the re¬ 
ceiver side, after a proper processing, two separate ES-MMSE 
equalizers and symbol-by-symbol receivers are adopted. Un¬ 
like the previous scenarios, in the presence of nonlinear distor¬ 
tions and phase noise the Alamouti scheme cannot perfectly 
separate the two signals at the receiver. However we will show 
in the numerical results that its performance is still excellent. 

The complexity of the channel model does not allow (to 
the best of our knowledge) to obtain results in a closed 
form as in previous sections. Hence, the achievable SEs for 
this scenario are computed through the Monte Carlo method 
proposed in |fT9) (see also pO) for details on the application 
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to the satellite channel). We point out that these values are a 
lower bound to the actual SE and they are achievable with the 
specific adopted receiver. The ensuing SE curves describe the 
possible achievable gains in a scenario that is more realistic 
than those described in previous sections. All results will be 
reported as function of Pgat /N, where Pgat is the HPA power 
at saturation. 


Modulation 

PaatjN [dB] 

QPSK 

-10 - 0 

8 PSK 

0 - 7.5 

16APSK/16PSK 

7.5 - 25 


Table I 

Psat/N RANGE OF THE ENVELOPE /j_p FOR PSK AND APSK 
MODULATIONS AND 7=1. 


A. Numerical results 

We consider transmitted signals with baudrate 37 Mbaud, 
adopting the classical constellations of satellite communica¬ 
tions, i.e., QPSK, 8PSK, 16APSK and 32APSK (denoted to 
as PSK/APSK schemes). As an alternative to the classical 
PSK/APSK, we also consider the use of 16PSK and 32PSK, 
as suggested by the theoretical analysis. The adopted shap¬ 
ing pulse p{t) has root raised cosine spectrum with roll-off 
a = 0.1. The input back-off is set to 0 dB for QPSK and 
8PSK, and to 3 dB for all other modulations]^ 


Pig. 15 shows the envelope of the pragmatic SE Ij p for 


the considered modulations, with power unbalance 7^ = 0 
dB. Details on the modulations of the envelope are reported 
in Table [Ij The figure also shows the SE for EDM, for the 
Alamouti scheme, and for a single satellite. In case of EDM, 
each signal has baudrate 1 /T =18.5 Mbaud, and the frequency 
spacing is equal to /c = (1 + ct)/T= 20.35 MHzj^ We can 
see from the figure that two overlapped signals can achieve 
a higher SE than all considered alternatives. Moreover the 
envelopes show that 32APSK and 32PSK are not convenient 
in case of overlapped signals, since they perform worse than 
16APSK and 16PSK modulations, and PSK modulations per¬ 
form better than APSK modulations. It is interesting to notice 
that, although the channel model is affected by nonlinear 
effects, inequality 0 still holds true even in this case. We 
also notice that, at high Pgat/N, EDM performs worse even 
than a single satellite. This loss is due to the interchannel 
interference (ICI) from the second EDM signal, which lies 
in the same OMUX bandwidth. In fact, due to the spectral 
regrowth after the HPA, the two EDM signals are no more 
orthogonal. This effect is proved in Pig. that compares the 
EDM curve with two SE curves: ideal EDM in the absence of 
ICI, and a single satellite with twice the power Psat- Similarly 
to the linear channel, ideal EDM can achieve the same SE as 
the single satellite with double power, but in the actual case 
ICI has an impact on performance. 


We can notice from Fig. 15 that gains given by two over¬ 
lapped signals w.r.t. a single satellite can be higher than 3 dB. 
The gains over 3 dB are related to the shaping of the overall 
signal, obtained by the sum of the satellite outputs. Indeed, as 


already mentioned in Section IV the sum of two signals has an 


amplitude distribution that is closer to a Gaussian distribution. 


Figs. 17 and 18 show the PDF and the cumulative distribution 
function (CDF) of the signal amplitude, properly normalized 


^We found these values to be optimal from other activities beyond this 
paper. We also point out that the impact of interchannel interference due 
to transponders transmitting on adjacent frequencies is negligible for all the 
presented scenarios, and hence it will not be considered | 21 |. 

^Other values of frequency spacing have been tested, but 20.35 MHz has 
been found to be practically optimal for this scenario. 



Figure 15. Spectral efficiency achievable by PSK and APSK constellations 
for two satellites and 7 ^ = 0 dB. 


by the number of transmitting satellites. We compare the 
amplitude distributions of a single signal and two overlapped 
signals (with 7^ = 0 dB), when the transmitters adopt 16PSK, 
RRC pulses with roll-off a = 0.1, IBO equal to 3 dB. For 
comparison purpose we report also the PDF and CDF of the 
Gaussian distribution with unit variance. It is clear from the 
figures that the sum of two signals is closer to a Gaussian 
distribution than the single transmitter. We have verified that 
similar considerations hold for 8PSK. 

Fig. 19 and Table [n| report SE curves for the same scenario 
as in Fig. 15 but with power unbalance equal to 6 dB. 



Figure 16. Spectral efficiency achievable by PSK/APSK constellations for 
two satellites using FDM and 7 ^ = 0 dB. 
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Signal amplitude 



P,at/iV / (dB) 


Figure 17. Probability density function of the signal amplitude with 16PSK, Figure 19. Spectral efficiency achievable by DVB-S2 constellations for two 
RRC a = 0.1, IBO=3 dB, 7 ^ = 0 dB. satellites and 7 ^ = —6 dB. 



Signal amplitude 


Figure 18. Cumulative distribution function of the signal amplitude with 
16PSK, RRC a = 0.1, 1B0=3 dB, 7 ^ = 0 dB. 






i?i / (bit/s/Hz) 


Figure 20. Spectral efficiency regions for DVB-S2 constellations for two 
satellites and 7 ^ = — 6 dB. 


Overlapped signals again outperform FDM for every Pgat/^ 
value but, since the equal rate constraint limits the performance 
to that of the lower power signal, the Alamouti scheme and a 
single satellite have higher SE at low Psat/N. The behavior 
of /j p w.r.t. a single satellite can be seen from the SE regions 
in Fig. 1^ and it can be noticed that it is perfectly in line with 
results found for the peak limited AWGN channel, despite a 
huge difference between the two models. 


Modulation 

Psat/N [dB] 

QPSK 

-10 - 5 

8 PSK 

5 - 12.5 

16APSK/16PSK 

12.5 - 25 


Table II 

Psat/N RANGE OF THE ENVELOPE /j,p FOR PSK AND APSK 
MODULATIONS AND 7=1/2. 


IX. Conclusions 

We investigated the rates achievable by a system using two 
co-located satellites. We exploited the second satellite to im¬ 
prove the spectral efficiency. We studied three models: AWGN 
channel with average power constraint, AWGN channel with 
peak power constraint, and the DVB-S2 satellite channel. For 
all cases we considered signals with overlapping frequencies, 
FDM, and the Alamouti scheme. Overlapped signals resulted 
to be convenient in all cases w.r.t. FDM, but we showed that 
there are cases in which the Alamouti scheme can outperform 
both, and that even a single satellite can be convenient over 
overlapped signals: these cases depend on the power unbalance 
and on the received signal-to-noise power ratio. 

Appendix A 

Alamouti Scheme with Time Misalignment 

In this appendix, we show an alternative implementation of 
the Alamouti scheme, for the case when the signals received 
by the two satellites have a time misalignment. The precoding 
and decoding schemes adopt complex conjugation and time 
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reversal of the signals. A similar precoding has been pro¬ 
posed in for orthogonal frequency division multiplexing 
schemes with two relay nodes. 

Let us consider a channel for which, when Xi{t) and X 2 {t) 
of hnite duration are transmitted, the received signal is 


obtained by sampling the continuous waveforms. The observ¬ 
ables for the two subchannels are 

yi = xi+ wi (23) 

2/2 = 72^2 + W2 (24) 


ya{t) = xi{t) + ^e^'^X 2 {t -t) + Wa{t) 

where r denotes the time misalignment. After the transmission 
of these signals, the transmitter sends x^i—t) and —x\{—t). 
The received signal, in this case, is 

yb{t) = x^i-t) - + t) + Wh{t) . 


If the receiver has perfect knowledge of r, (p, and 7 , it can 
elaborate the received signals as 


Vait) - + t) 

\/l + 7^ 
v/i + 7^ 


\/l+ 'y'^Xi{t) + Wait) 
yrr l‘^X2it) + Wbit) , 


where the Gaussian processes Wait) and Wbit) are statistically 
equivalent to Wait) and Wbit). Hence, signals xiit) and X 2 (f) 
can be independently detected. 

In the presence of phase noise (i.e., when the phase p is not 
constant) and nonlinear effects, however, this scheme works 
only approximately. The loss due to residual interference will 
be negligible if signals xi(f) and X 2 (f) have a duration shorter 
than the phase noise coherence time and nonlinear effects 
are limited. In particular, the AM/PM characteristic of the 
nonlinear ampliher must be such that when conjugating the 
input, the output results to be conjugated too. 


Appendix B 
Proof of Theorem[T] 

We now prove Theorem [T] we first prove a preliminary 
result concerning the differential entropies of two continuous 
random variables. 

Lemma. Let x and y be two independent continuous complex 
random variables, with probability density functions pix) and 
piy) and differential entropies hix) and hiy). Then 

+ + ( 20 ) 

with equality if and only if x and y are independent Gaussian 
random variables with the same variance. 


Proof: For the entropy power inequality 

2h{x+v) ^ 2 ^(^) _|_ 


= 2 


1 - 1 - 


h{x)-\-h{y) 


cosh 


hix) - hiy) 


In 2 


> 2 


i-t 


h(x) + h(y) 


( 21 ) 


( 22 ) 


where equalities in and ( | 2 ^ hold if and only if x and 
y are Gaussian and have same variance. Eq. ( | 20 l i is hnally 
derived by taking the logarithm of ( | 2 ^ . ■ 

We then consider the rates achievable by EDM. Under the 
assumption of ideal EDM transmission, a sufficient statistic is 


where xi and X 2 are the signal samples, rci and W 2 are white 
Gaussian noise processes with power N/2 instead of N, since 
EDM works with half the bandwidth w.r.t. the case of a single 
transmitter. The mutual information of EDM is the average of 
the mutual information for the two channels, i.e.. 


IpDM 


hjyi) + ^( 2 / 2 ) 
2 


l 0 g 2 



and the pragmatic rate is 

Ffdm.p = hiy 2 ) — log2 

Since the mutual information is a non decreasing function of 
the SNR 1^, clearly it is /fdm,p < Ffdm- 
We finally prove Theorem [T] 

Proof: We first prove inequalities (|^ and 



2/2 > /fdm,p ■ (25) 

The samples at the output of channel Q can be equivalently 
expressed as ?/ = 2 /i + 2/2 the mutual information of this 
equivalent expression reads /j = /i(j/i -f 2 / 2 ) — log 2 ( 7 reAf). 
Hence, 

from which, by an application of the Lemma, we derive 
inequality Q. The mutual information I 2 , instead, reads 

I 2 = hiy\xi)-log^iTteN) 

= hiy 2 Pwi)- log 2 iireN) 

and 


2/2 - f^FDM,p = 2 /i(?/2 + wi) - h(2/2) - log2(27reA^) 

which becomes ( |25] l from Lemma. 

Since /j_p = min (/j, 2 I 2 ) and /fdm > Ffdm.p, clearly 
follows with equality if and only if xi and X 2 are Gaussian 
with 7 ^ = 0 dB. ■ 


Appendix C 
Proof of Theorem|2] 

We now prove Theorem]^ 

Proof: Let us start by hrst proving inequality (a). The 
observable for Alamouti precoding is 


2/A,i = -I- 7 x 2 -f WA,i 

2/A.2 = X2- yxi -I- WA.2 

where wa,i,wa ,2 are independent Gaussian random variables 
with power N. 
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Let us evaluate 



1 

II 

1 

h{yA.i,yA,2) 

(26) 

2 

IV 

1 

HyA,i) + h{yl2) 

2 

(27) 

= 0 


(28) 


where ( |28] l is obtained by observing that h{y) = h{yx,i) = 
h{y\ 2 )- Equality in ( |27| l is achieved if and only if ?/a,i and 
y\ 2 are independent. From an application of Lukacs-King 
theorem p4| , independence holds if and only if xi and X2 are 
independent Gaussian random variables with the same power. 

We now prove inequality (&); the Alamouti observable, after 
receiver processing, is 

yA,i = x/TT^Xi + r&A.i i = 1,2. (29) 


with SNR alP/N. Say I{alP/N) the mutual information 
of (|3 T]i as function of the SNR, the information rate is 



where the last inequality is due to the concavity of the mutual 
information as a function of the SNR j^ . The trace can be 
rewritten as 

trace h'^ = trace (l + 2 -fn{^H) + . 

Since we are dealing with an ergodic process, it holds that 


and is still a sufficient statistic for detection. The SNR in ( |29| ) 
is ( 1 + 7^)/2 < 1 times the one in ( | 2 ^ and (l + 7 ^)/ 27 ^ > 1 
times that in Hence, since the mutual information is 

a concave function of the SNR no matter the input 

distribution, inequality (b) is straightforward and it holds with 
equality if and only if 7 ^ = 0 dB. ■ 

Appendix D 

Ergodic Rates for Satellites Transmitting the 

Same Signal 

We demonstrate that, when the two satellites transmit the 
same signal, the information rate cannot be higher than the 
case when transmitting independent signals. 

Proof: When the same signal is transmitted, we obtain 
a time-varying frequency selective channel as in (|7]). The 
received signal can be expressed, by using a matrix notation, 
as 

y = Hx + w , (30) 

where x = [xq, ..., is the vector of transmitted 

symbols, w is the noise vector, and H is a matrix with 
elements Hk^i = In this section we denote the Hermitian 
operator by f, and the identity matrix by I. 

Matrix H, for the channel we are considering, can be 
written in the following form 

H = I + -f^H, 

where matrix $ is diagonal with elements €>fc,fc = and 
H is Toeplitz with elements Hk,i = sinc((j — k)T — r). 

We denote the singular value decomposition of H as 

H = , 

where U, are unitary matrices, and S is a diagonal matrix 
with elements 'Ek,k = crk- 

If we set X — X and y = U^y, channel ( [30l l becomes n 
parallel channels in the form 

yk = cFkXk+Wk, (31) 


lim —trace {h^h\ = lim —trace 

n—)-c 50 n \ J n—>-oo Ji \ 

= 1 + 7^ 


where last equality is found through the Szego Theorem | [25j . 
We finally obtain the following inequality 


lim 

n—foo 


n 


< 


I 



The right term is independent of U and . Hence, two satel¬ 
lites transmitting the same signal have achievable information 
rate lower than the one achievable with the Alamouti scheme 
and same input distribution. ■ 
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